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ABSTRACT 

We present Far Ultraviolet Spectroscopic Explorer (FUSE) and Space Tele- 
scope Imaging Spectrograph observations of the weak interstellar N I A1160 dou- 
blet toward 17 high-density sightlines [N{}itot) > 10^^ cm~^]. When combined 
with published data, our results reveal variations in the fractional N I abundance 
showing a systematic deficiency at large N(B.tot)- At the FUSE resolution (~ 
20 km s~^) the effects of unresolved saturation cannot be conclusively ruled out, 
although O I A1356 shows little evidence for saturation. We investigated the 
possibility that the N I variability is due to the formation of N2 in our mostly 
dense regions. The 0-0 band of the c'^^^u ~ -^^^t transition of N2 at 958 A 
should be easily detected in our FUSE data; for 10 of the denser sightlines N2 
is not observed at a sensitivity level of a few times 10^^ cm~^. The observed N I 
variations are suggestive of an incomplete understanding of nitrogen chemistry. 

Subject headings: ISM: atoms — ISM: abundances — ISM: clouds — ultraviolet: 
ISM 

1. Introduction 



Elemental abundance studies are important for models of Galactic chemical evolution. 
Studies of the interstellar oxygen abundance support a relatively constant value of O/Rtot = 
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(3.43 ± 0.15) X 10 ^ out to ~ 1000 kpc, with any variability less than the 1 a measurement 
uncertainties (Meyer, Jura, & CardeUi 1998; Cartledge et al. 2001; Andre et al. 2003). The 
situation is not as clear for carbon and nitrogen. In this work, we focus specifically on the 
abundance of interstellar nitrogen. 

It is generally believed that interstellar nitrogen is a product of the CNO cycle and 
recycled into the ISM through the winds of low and intermediate mass stars (e.g., red giant 
branch and/or asymptotic giant branch stars (AGB); Pilyugin, Thuan, & Vichcz 2003). In 
contrast, oxygen is produced in stars during the He burning phase and is returned to the 
ISM via supernovae. 

Early studies (Lugger et al. 1978; Fcrlct 1981; York et al. 1983) of interstellar N I 
using the Copernicus satellite investigated moderately dense lines of sight and found that 
A^(N I) increases linearly with N{Y[tot) [N{B.tot) = 2A^(H2) + iV(H I)]. Using high-quality 
Goddard High-Resolution Spectrograph (GHRS) data for seven sight lines with Nljitot) < 
10^1 cm ^, Meyer, CardeUi, & Sofia (1997) suggested that the interstellar nitrogen abundance 
is constant, N/H^t = (7.5 ± 0.4) x 10"^. However, Jenkins et al. (1999) and Sonneborn 
et al. (2000), using data from the Interstellar Medium Absorption Profile Spectrograph 
(IMAPS), show a factor-of-2 variation in the N I abundance between 5 Ori [N l/Rtot — (3.97 
± 0.30) X 10-5] and 7^ Vel [N l/Rtot = (7.99 ± 0.47) x 10"^]. 

In the near interstellar medium (ISM), d < lOOpc, N I is more affected by ionization 
than O I because O I is coupled more strongly to H I by charge exchange reactions (Sofia 
& Jenkins 1998; Jenkins et al. 2000; Lehner et al. 2003). Therefore, these near ISM 
sightlines (Lehner et al. 2003) will not be discussed here. The effects of ionization at larger 
column densities are thought to be insignificant. In addition, nitrogen is not incorporated 
into refractory interstellar dust grains (Sofia, CardeUi, & Savage 1994). Nitrogen bearing 
interstellar ices (Gibb, Whittet, & Chiar 2001; Chiar ct al. 2002) arc not expected along 
our moderately reddened {Ay < 2.0) sight lines since interstellar ices do not form until A^, 
> 3 (Whittet et al. 2001). However, nitrogen chemistry and N2 formation should become 
important. 

In order to probe the extent of potential variations in the interstellar nitrogen abun- 
dance, we have undertaken a survey of the weak interstellar N I doublet AA1159.817, 1160.937. 
Use of this doublet removes uncertainty in oscillator strengths (/-values) for different tran- 
sitions and allows for a direct comparison to the work of Meyer et al. (1997). Our survey 
utilized archival data from the Far Ultraviolet Spectroscopic Explorer (FUSE; Moos et al. 
2000j and the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space 
Telescope (HST) toward 17 high-density [A^(H) > 10^^ cm sightlines. 
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2. Observations and Data Reduction 
2.1. FUSE Data 

The N I data for 16 of the stars studied here, shown in Table 1, were obtained throughout 
the FUSE Prime mission phase (November 1999 — March 2003). All of the data were 
acquired with the star in the large aperture (30" x 30"; LWRS) with the exception of 
HD 219188 for which the medium aperture (4" x 20"; MDRS) was utilized. The data cover 
the wavelength range 905 — 1185 A with a spectral resolution equivalent to At; ~ 20 km s~^. 
The weak interstellar N I doublet at 1160 A appears in both the LiFl and LiF2 channels 
which provides a consistency check to rule out the possibility of detector artifacts. 

The time-tagged and histogram data were reduced and calibrated with CalFUSE^ (ver- 
sion 2.2.2; Dixon & Sahnow 2003). CalFUSE provides the appropriate Doppler corrections 
to remove the effects of the spacecraft motion and places the data on the heliocentric ve- 
locity scale iyheiio)- The wavelength solution provides good relative calibration across the 
LiF channels. In order to minimize the uncertainties in the relative wavelength calibration 
between exposures, the data were co-added with a cross correlation technique. We used in- 
terstellar H2 hues in several of the FUSE channels to fix the velocity scale (KeO- The final 
summed spectra (near 1160 A) have signal-to- noise (S/N) ratios between ~ 20 and ~ 300 
per resolution element for all data sets. S/N ratios greater than 30 were obtained utilizing 
focal plane splits or other similar procedures. 



2.2. STIS Data 

STIS observed HD 147888 (N I A1161) for 1656 s on 2000 Aug 17, HD 24534 (O I A1356) 
for 2945 s on 2001 Mar 25, and HD 219188 (O I A1356) for 1200 s on 2001 Sept 27. The N I 
A1161 data for HD 147888 used the E140H grating centered at 1271 A and the 0.2" x 0.09" 
aperture resulting in a resolving power {R) of about 100,000 (At; ~ 2.8 km s~^). The O I 
A1356 data utilized the same grating but with the 0.1" x 0.03" aperture centered on A1416 
for HD 24534 and on A1271 for HD 219188 resulting in - 200,000 (Av - 1.5 km s"^). 

The E140H setup centered at A1271 used to measure N I toward HD 147888, is not ideal 
since the stronger N I line is not obtained. In addition, the MgF2 detector windows severely 
attenuates the flux below 1200 A. The 1161 A line can only be detected in relatively dense 
sightlines or with long exposure times. With the exception of HD 147888, examination of 



•^The CalFUSE Pipeline Reference Guide is available at http://fuse.pha.jhu.edu/analysis/pipeline_reference.litnil. 
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the archival data for other sighthnes did not reveal evidence for the N I 1161 A line due to 
the relatively low S/N (< 15 per resolution element). 

The data were reduced and extracted with the CALSTIS pipehne (v2.13b). The sub- 
traction of background and scattered hght from the echelle data employs the algorithm of 
Lindler & Bowers (2000). CALSTIS provides the appropriate Doppler corrections to place 
the spectra on the Vheiio scale. The spectrum of HD 147888 at N I A1161 has a S/N ratio of 
~ 12 per resolution element. For O I A1356 toward HD 24534 and HD 219188, we obtained 
a S/N of ~ 70 and 48 per resolution element, respectively. 

3. Results and Discussion 

Figure 1 shows the apparent column density (ACD) profiles, Na{v), for both members of 
the weak interstellar N I doublet toward 2 stars of our stellar sample. See Savage & Sembach 
(1991) for a detailed discussion of the apparent optical depth method. The excellent overall 
agreement between the Na{v) profiles for both members of the N I doublet implies that the 
apparent column densities are the true column densities. Table 1 presents our N I column 
densities, the average of both members of the doublet except HD 147888 (assumed to be 
unsaturated). Also presented are the measured equivalent widths (Wxs) of the N I doublet 
for each star. Finally, Table 1 shows the available Copernicus, International Ultraviolet 
Explorer, GHRS, IMAPS, and STIS measurements of A^(N I), N(litot), and N{0 I). 

Figure 2 {Top) shows A^(N I)/A^(Htot) versus N(H.tot) which reveals a departure from the 
average at high column densities. With the exception of 5 Ori, the N I abundance appears to 
be constant for 7V(Htot) < 10^^ cm~^, but an anticorrelation seems to be present for greater 
column densities. A Student t test yields a 7.6% probability that the data sets for N{}itot) 
< 10^^ cm~^ and NCHtot) > 10^^ cm~^ (including the outlier 5 Ori) arise from the same 
parent population. If S Ori is excluded, the probability is reduced to 1.3%. Therefore, we 
conclude that the observed variability is real. Figure 2 {Bottom) depicts A^(N l)/A^(0 I) as 
a function of N{lltot)- The data point for HD 147888, the densest sightline presented here, 
is high because O I is slightly depleted (Cartledge et al. 2001). The large scatter in N/0 
supports that the variation in N I is not an artifact introduced by the larger uncertainties 
for N{}lt,t). 

In order to test whether the observed N I variability is due to observational effects (i.e., 
saturation) we calculated curves-of-growth for four stars which exhibited deviant behavior 
(e.g., HD 179406) and six stars which did not. The curve of growth results agree with those 
obtained from the ACD method. Additionally, our results compare favorably to previous 
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N I measurements (Hoopes et al. 2003; Sonnentrucker 2003, private communication). We 
believe saturation effects to be minimal since tn i/tq i = 1.4 and O I shows little evidence 
for saturation (e.g., Cartledge et al. 2001). 

A second source of systematic error could arise if the instrinsic absorption profile con- 
sisted of a narrow, but unsaturated, feature plus broad shallow wings. In such a scenario, 
the weak N I doublet associated with the broad component could be too weak to be detected 
at low S/N. If enough H I column is located in the broad component, this could skew the 
N l/Htot ratio to systematically low values. We can investigate this possibility in two ways: 
1.) through a comparison of the measured Doppler broadening parameter, 6- value, for the N I 
and hydrogen lines and 2.) comparing N 1/B.tot vs. distance to test whether systematically 
lower ratios are caused by additional weak kinematic components for the longer lines-of-sight 
(Spitzer 1985). 

Unfortunately, both the Ly-a line of H I and the J = 0, 1 lines of H2 typically lie 
on the square-root part of the curve of growth for these sightlines and are insensitive to 
the &- value. However, &- values can nominally be measured for the J > 2 lines of H2. The 
pubhshed 6-values for HD 73882, HD 110432, HD 185418, and HD 192639 (Rachford et 
al. 2001; Sonnentrucker et al. 2002; Sonnentrucker 2003, private communication), all show 
comparable or smaller 6-values for H2 than for N I indicating that all the N I is detected. 
Additionally, no statistically significant slope of N l/Htot vs. stellar distance (125pc < d < 
2kpc) is found for our sample. These points argue against intrinsic line shape contributing 
to the N 1/B.tot deficiencies reported herein. 

3.1. N2 

Based on models of steady-state, gas-phase interstellar chemistry, N2 is expected to be 
the most abundant nitrogen-bearing molecule in dense clouds. Viala (1986) predicts that 
at the column densities of the sight lines studied here [N{}itot) of a few times 10^^ cm^^], 
the N2 column densities should be on the order of 15% of the N I abundance or A^(N2) ~ 
10^^ cm~^. Inclusion of time dependence or depletion onto grain mantles (Bergin, Langcr, 
& Goldsmith 1995) could result in a slightly smaller gas-phase N2 abundance. However, N2 
is still predicted to be the most abundant nitrogen bearing molecule. 

The strongest band of N2, covered by FUSE, is the 0-0 band of the C4^E+ — X^S+ 
transition of N2 at 958 A. Other N2 bands reside in the far- ultraviolet, but are blended 
with more abundant species (e.g., H2). Utilizing the laboratory wavelengths and /-values 
for the 0-0 band (Stark et al. 2000), we created several synthetic N2 spectra, assuming level 



-6- 



populations for excitation temperatures between 10 — 1000 K (c.f., McCandliss 2003). Our 
synthetic spectrum is saturated for N(N2) — 10^^ cm~^. Therefore, we searched our FUSE 
data for the presence of N2. 

Figure 3 (Top) exhibits the 958 A portion of the spectrum toward HD 210839 [E{B 
— V) — 0.62 mag.] and the bottom panel shows our synthetic N2 spectrum for A'"(N2) = 
10^^ cm~^, 10% of that predicted. The lower column density incorporates the effects of four 
velocity components detected in O I A 1356 (Andre et al. 2003) along the line of sight with 
similar amounts of material. There is no evidence for N2 in the spectrum of HD 210839 to a 
level of a few times 10^^ cm~^. Examination of the ten densest sightlines that have FUSE 
data also show a similar dearth of N2. Hence, gas-phase N2 cannot explain the observed N I 
variability. Quantitative limits on N2 will be presented in a subsequent paper. 

4. Summciry 

Our data show that the interstellar N I abundances relative to H I and O I appear to vary. 
Very high-resolution studies of these and other sight lines will help to verify and expand the 
findings reported here. Although N2 is expected to be the most abundant nitrogen-bearing 
molecule in moderately dense [log A^(H2) > 10^^ cm~^] clouds (e.g., Viala 1986), our data do 
not support the predicted abundances. The presence of an anticorrelation above N(lifot) — 
10^^ cm~^ suggests that the mostly likely explanation of the interstellar N I variability is the 
need for a better understanding of nitrogen chemistry. However, differences in the mixing 
processes of products from AGB stellar outflows and SN may also contribute. 

We thank the FUSE Science Operations team for their dedication and efforts in acquiring 
these observations. This research made use of the Simbad database, operated at CDS, 
Strasbourg, France. 
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Fig. 1.— Top: N I apparent column density profiles for HD 191877 and HD 195965. The 
stronger line is denoted by the solid circles, while the weaker line is represented by the open 
squares. The vertical lines represent 1 a uncertainties that were determined following the 
procedures outlined in Sembach & Savage (1992). The signal-to-noise ratios are 320 per 
resolution element for HD 191877 and 270 per resolution element for HD 195965. There is 
no evidence for unresolved saturated structure. 
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Fig. 2.— Top: N(N l)/N{B. 

tot) cis a function of N[}ltot)- The dashed hne represents the 
average value N/Htot from Meyer et al. (1997). For N{B.tot) > 10^^ cm~^ there is a large 
dispersion from the average, unlike that found for interstellar O I (e.g., Meyer et al. 1998). 
The observed variability either represents a real cosmic variance or evidence for enhanced 
chemical processing of N I. Bottom: N(N l)/N{0 I) as a function of N{B.tot)- The dashed 
line is Meyer's average values for N/0 [i.e., (N /Rtot) / {O /Rtot)', Meyer et al. 1997; Meyer et 
al. 1998]. The large scatter verifies that N I varies with N{B.tot) and is not an artifact of the 
larger uncertainties associated with N(lltot)- 
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Fig. 3.— Top - The 958 A portion of the FUSE data toward HD 210839 (S/N ~ 30 per 



resolution element). Bottom — A synthetic spectrum of the 0-0 band of the C4 



transition of N2 at 958 A, which is the strongest N2 band in the FUSE wavelength region. 
The column density of the synthetic spectrum is N{N2) — 10^^ cm~^, 10% of that expected 
from steady-state gas- phase chemistry (Viala 1986). The synthetic spectrum was produced 
at T = 100 K, b = 5 km (narrow histograms) and convolved to the FUSE instrumental 
resolution ~ 17,000 (solid line). N2 should be easily detected in the FUSE data but it is 
not. 
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Table 1. Interstellar N I, O I, and N I/O I Ratios^ 
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5 














11, 15 


HD 94493 


7.7 ± 1.2 


2.1 


± 


0.7 


(7.64 


± 


1.93) 


X 


1016 


(4.52 ± 


1.21) 


X 


10- 


-5 


(7.0 


± 


0.6) X 10" 


0.11 


± 


0.03 


11, 13 


HD 24534 


10.2 ± 0.9 


3.1 


± 


0.8 


(1.10 


± 


0.20) 


X 


10" 


(5.02 ± 


1.25) 


X 


10- 


-5 


(7.43 


± 


0.28) X 10" 


0.15 


± 


0.03 


11, 15 


HD 185418 


16.0 ± 2.0 


2.7 


± 


1.5 


(1.35 


± 


0.37) 


X 


10" 


(5.42 1 


1.94-, 
1.73/ 


X 


10- 


-5 


(1.20 


+ 


°0;f2) X 1018 


(0.11 


+ 


0.03 
0.03'' 


11, 13, 14, 16 


HD 99857 


15.5 ± 1.3 


5.2 


± 


0.9 


(1.78 


± 


0.24) 


X 


10" 


(7.09 ± 


1.14) 


X 


10- 


-6 


(7.8 


± 


0.5) X 10" 


0.23 


± 


0.03 


11, 13 


HD 210839 


16.5 ± 0.7 


5.4 


± 


0.8 


(1.84 


± 


0.19) 


X 


10" 


(6.53 ± 


1.00) 


X 


10- 


-5 


(1.30 


± 


0.06) X lOis 


0.14 


± 


0.02 


11, 13 


HD 179406 


10.5 ± 1.4 


3.8 


± 


1.3 


(1.25 


± 


0.31) 


X 


10" 


(4.15 ± 


1.01) 


X 


10- 


-5 














8, 9, 11 


HD 192639 


12.5 ± 2.7 


3.3 


± 


1.6 


(1.30 


± 


0.46) 


X 


10" 


(4.31 ± 


1.64) 


X 


10- 


-T, 


(1.35 


± 


0.06) X lOi* 


0.10 




0.03 


11, 13, 17 


HD 124314 


19.5 ± 1.2 


6.5 


± 


1.1 


(2.21 


± 


0.27) 


X 


10" 


(6.95 ± 


1.21) 


X 


10- 


-5 


(1.53 


± 


0.08) X lOis 


0.14 


± 


0.02 


11, 13 


HD 73882 


16.0 ± 1.6 


5.6 


± 


1.6 


(1.88 


± 


0.38) 


X 


10" 


(4.87 t 


2.54-, 
1.86/ 


X 


10- 


-5 














11, 15 


HD 147888 




7.3 


± 


1.5 


(3.41 


± 


0.78) 


X 


10" 


(5.78 ± 


1.59) 


X 


10- 


-6 


(1.41 


± 


0.11) X 10i8 


0.24 


± 


0.06 


11, 14 
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